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Simple sets of equations have been 
derived to describe the absorption of 
microwave power in three-region, lossy 
transmission lines in terms of S-parameter 
reflection and transmission ampHtudes. 
Each region was assumed to be homo- 
geneous with discontinuities at the region 
boundaries. Different sets of equations 
were derived to describe different 
assumptions about the amplitudes of the 
reflection coefficients at the different 
boundaries. These equations, which are 
usefiil when interference effects due to 
multiple reflections are small, were used to 
analyze S-parameter measurements on a 
transmission line that had a microfluidic 
channel in its middle region. The channel 
was empty for one set of measurements 
and filled with water for a second set of 
measurements. Most of the reflection 
assumptions considered here produced 
similar results for the traction of the 
applied microwave power that was 
absorbed by a water-filled microchannel. 
This shows that the absorbed power is 
relatively insensitive to the reflection 
details as long as energy is conserved in 



the analysis. Another important result of 
this work is that the difference between 
the power absorbed in a water-filled 
channel and the power absorbed in the 
same empty channel can be a poor 
predictor of the power absorbed in the 
water in the presence of competing 
absorption processes such as absorption 
by the transmission-line metal. 
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1. Introduction 

The possibility of microwave heating of lossy fluids 
in microfluidic devices with integrated transmission 
lines has been demonstrated recently [1, 2, 3, 4]. To 
optimize the performance of microwave heating in 
these types of devices, it will be necessary to distin- 
guish between the absorption of microwave power in 
fluids in the microchannel and the absorption of 
microwave power in the metal components of the trans- 
mission line. 

The standard way to calculate the power absorbed in 
a transmission line [5] requires knowledge of the real 
and imaginary components of the propagation con- 



stants of the transmission lines connected to the trans- 
mission line being modeled. But it is possible that this 
type of information will not be available for transmis- 
sion lines that are integrated with microfluidic devices. 
Therefore, it may be desirable to be able to analyze 
S-parameter measurements of this type of device with- 
out access to the type of data needed for the 
model described in [5]. 

In general, the reflection and transmittance spectra of 
a complex transmission line will have interference 
fringes due to inter-reflections between discontinuities 
in the transmission line. However, the fringes will be 
small perturbations on the spectra when all of the 
reflection coefficients are small, which is usually a 
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design goal. In this case, it should be possible to ignore 
the phase information necessary to describe the ampli- 
tude and spacing of the interference fringes. 
The purpose of this paper is to develop approximations 
for the microwave power absorbed in lossy components 
of a three-region transmission line from S-parameter 
reflectance and transmittance amplitude measurements. 
The results should prove useful when the amplitude of 
the fringes on the absorption spectra are small com- 
pared to the effects of interest, as illustrated by applica- 
tion to experimental data on a real device. 

Section 2 of this paper derives equations for the frac- 
tions of power reflected from, transmitted through, and 
absorbed in each region of a transmission line made up 
of three different homogeneous regions with a single 
reflective boundary. In other words, the reflection coef- 
ficients at all but one of the region boundaries are 
assumed to be zero. Section 3 derives equations for the 
same quantities when symmetric pairs of the four 
boundaries are equally reflective. These equations 
express the S-parameter data in terms of attenuation 
coefficients within the different regions of the transmis- 
sion line and the reflection coefficients at the bound- 
aries of the transmission-line regions. Section 4 inverts 
these equations to express the reflection and attenuation 
coefficients in terms of the S-parameter data. The 
resulting equations are then used to analyze experimen- 
tal data from a three-region transmission line that 
contains a microfluidic channel in its middle region. 
Section 5 discusses the results and draws some con- 
clusions. 



2. Absorption in Transmission Lines 
With a Single Reflective Boundary 

Figure 1 illustrates a complex generic transmission 
line that consists of three regions of length Zj, z-^, and z^, 
respectively, separated by boundaries labeled A, B, C, 
and D. The geometry and dielectric constants of all of 
the structures that make up each of the regions are 
assumed to be uniform within each region. It is also 
assumed that microwave power propagates along the 
transmission line with no loss due to emission of 
radiation. Finally, it is assumed that Regions 1 and 3 
have identical microwave properties and the two trans- 
mission lines (not shown) to which the microwave 
probes are connected have identical microwave 
properties. 

Let P be the power incident on Boundary A from the 
left, and let RP be the power leaving Boundary A to the 
left. Similarly, let TP be the power transmitted to the 
right from Boundary D. Now assume that R and T, 
which will be called the reflection and transmission, 
respectively, are the only experimental data available 
for the device of Fig. 1 . Further assume that the real and 
imaginary parts of the propagation constants of the 
transmission lines (not shown) connected to the device 
shown in Fig. 1 , as well as those of the three regions in 
that device, are not known. In this case, there is not 
enough information available to invert Maxwell's 
equations to determine the propagation and attenuation 
constant in the regions of the device of Fig. 1 and in the 
transmission lines to which it is attached. This remains 
true even if Regions 1 through 3 are assumed to have 
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Fig. 1. A complex generic transmission line consisting of three regions showing the power propagating in each direction at certain locations in 
the transmission line. 
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identical properties and the transmission lines attached 
to the device in Fig. 1 are assumed to be identical. 

In this case it is convenient to treat the transmission 
line in Fig. 1 as a power- flow problem [6]. A fraction 
RP of the incident power is reflected at Boundary A, a 
fraction AiP is absorbed in Region 1, a fraction T12P is 
transmitted across Boundary B, and a fraction R12P 
is also reflected at Boundary B. Similar descriptions 
apply to Regions 2 and 3, and, finally, a fraction TP 
is transmitted out the right side of the transmission 
line at Boundary D. By conservation of energy, 
A=A-^+ A2 + At,= \—R — T'k the fraction of the power 
incident on the transmission line that is absorbed by 
the transmission line. R and T can be obtained directly 
from S-parameter measurements, and A can be 
calculated from R and T. 

On the other hand, ^1, A2, and A-^ are complex func- 
tions of the microwave power-reflection coefficients 
and the phase differences between the microwave 
amplitudes reflected at and transmitted across the 
different boundaries. In general, these quantities are 
described by Maxwell's equations. However, if inter- 
ference between multiply reflected waves in the trans- 
mission line introduces only small effects (fringes on 
the absorption spectra), compared to the effects of 
interest, then simpler descriptions are possible. These 
are based on ray tracing, knowledge of the general 
nature of the solutions of Maxwell's equations in 
homogeneous lossy regions, conservation of energy, 
and a number of independent reflection and transmis- 
sion measurements equal to or greater than the number 
of unknown parameters. This idea is developed in the 
following two sections. 

2.1 Reflection Only at Boundary A 

First assume that the power-reflection coefficients at 
Boundaries B, C, and D are all zero. Then i?i2 = ^23 = 
i?3 = 0. In this case. 



r=[l-p]TiT2T3=10<^^l''°\ 


(1) 


R =p= 10^10), 


(2) 


A, = {\-p\{\-T,l 


(3) 


^2 = [1-P]T:[1-TJ, 


(4) 


A, = {\-p\XiX2{\-x,l 


(5) 


T„ =exp(-2a^„). 


(6) 



In these equations, Sn and S21 are the respective reflec- 
tion and transmission coefficients in dB, p is the power- 
reflectance coefficient S^ at Boundary A, and a„ for 
« = 1, . . ., 3 are the amplitude-attenuation coefficients 
for Regions 1-3, respectively. Even though the order of 
the factors in these equations is not important, it simpli- 
fies tracing the flow of power through the device. 

The amplitude absorption coefficient of a material, 
which is often denoted by a, should not be confused 
with the amplitude attenuation coefficient a„ of Region 
n of the transmission lines being considered here. Two 
times the amplitude absorption coefficient of a materi- 
al, which can be calculated from the complex dielectric 
constant of that material, describes the power attenua- 
tion per unit distance of a collimated wave of electro- 
magnetic radiation propagating through a homo- 
geneous slab of that material. On the other hand, two 
times a„ describes the power attenuation per unit 
distance of a guided wave of radiation propagating in 
Region « of a complex transmission line. In general a„ 
will depend upon both the geometry and the 
frequency-dependent complex dielectric constants of 
all of the structures that make up Region n of the trans- 
mission line. 

2.2 Reflection Only at Boundary B 

Only when all but one of the reflection coefficients 
are zero are simple equations like Eqs. (1-5) rigorously 
true without regard to the phase relations of the 
microwave amplitudes. For instance, let p be the power 
reflection coefficient at Boundary B, and let all of the 
other reflection coefficients be zero. In this case, 
Eq. (2) and Eq. (3) must be replaced by 

R =TiPTi, (7) 

A, = {\- T,] + T,p [1 - T,], (8) 

while the other equations remain valid. 

2.3 Reflection Only at Boundary C 

If p is the power reflection coefficient at Boundary C 
and all of the other reflection coefficients are zero, then 
Eqs. (2-4) must be replaced by 

R =TiT2pT2Ti, (9) 

^1 = [1 - Ti] + T1T2PT2 [1 - Tj , (10) 

A2 = Ti [1 - T2] + T1T2P [1 - T2] . (11) 
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2.4 Reflection Only at Boundary D 

Continuing in the same vein, if p is the power reflec- 
tion coefficient at Boundary D, and all of the other 
reflection coefficients are zero, then Eq. (1) is still 
valid, but Eqs. (2-5) must be replaced by 

R = T1T2T3PT3T2T1, (12) 

^1 = [1 - Ti] + T1T2T3PT3T2 [1 - Ti] , (13) 

A = Ti [1 - T2] + T1T2T3PT3 [1 - T2] , (14) 

^ = [1 - Ti] + T1T2T3P [1 - T3] . (15) 

3. Absorption in Transmission Lines 
Witli Two Reflective Boundaries 

In the case where the amplitude of any two of the 
reflection coefficients is non-zero, it is still possible to 
derive good approximations to the reflectance and 
transmittance based on power flow, if the effects of 
interference are small compared to the effects of inter- 
est. Practically speaking, this means that the reflectance 
coefficients will be small compared to one, which is 
equivalent to requiring good impedance matching at all 
boundaries. When the impedance matching is poor, 
there will be large fringes on the Sn and ^21 spectra due 
to interference between multiply reflected waves prop- 
agating on the transmission line. On the other hand, 
when the impedance matching is good, the fringes will 
be barely visible compared to the effects of interest in 
the spectra. 

In the approach used here [6], the powers in multiply 
reflected waves are simply added, which washes out all 
interference effects and produces a result intermediate 
between the envelope of the fringe maxima (construc- 
tive interference) and the envelope of the fringe mini- 
ma (destructive interference) in the S-parameter meas- 
urements. If the interference effects are large compared 
to the effects of interest, then the errors introduced by 
this approximation will be large. On the other hand, 
when the fringes are small compared to the effects of 
interest, then the derived equations will provide useful 
approximations. 

3.1 Reflection Only at Boundaries A and D 

Continue to assume that Regions 1 and 3 have iden- 
tical properties and that p is the power reflection coef- 
ficient at each probe. To derive a general expression, it 
is necessary to have expressions for the reflectance 
of the transmission line at the left probe and trans- 



mittance of the transmission line at the right probe. 
These can be derived [6] by expressing A^ A2, and A^ as 
well as R and Tas an infinite series of terms; each term 
describing the power in a wave that has been multiply 
reflected a different number of times, and summing the 
series analytically, to produce closed-form expressions: 
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^p[l+T^], 



1-p T 



[l-p][l-Ti][l+TiT2T3PT3T2] 
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-[l-Ti][l-p+TiT2T3PT3T2], 
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A = 



4 = 



(16) 
(17) 

(18) 

(19) 

(20) 
(21) 



The approximations on the right hand side of Eq. (16) 
and Eq. (20) were chosen to conserve energy as well as 
be accurate to order p^. They will he useful when 
p ^ « 1 . They could also have been derived by ray 
tracing [6] to generate all terms linear in p without 
recourse to the sum of the series of all ray-tracing 
terms. 

3.2 Reflection Only at Boundaries B and C 

As another alternative reflection model, assume 
again that Regions 1 and 3 have identical properties, 
but p is the power reflection coefficient at Boundaries 
B and C of Region 2. In this case, the only expressions 
that need to be derived by the procedure described above 
are those for Ri2 and ^23 , because Ai, A2, and A^ can be 
expressed in terms of these quantities as shown below: 
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p[l + T2^(l-2p)] 
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^p[l+T2'], (22) 
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7?=Ti7?j,T, »<p[l+T,'], (24) 

r = T7,3T3»TjT,T3[l-2p], (25) 

A,= [l-Tj[l+T,R,,], (26) 

A,=T,[l-R,,+T,^], (27) 

4=t7,3[1-T3]. (28) 

The approximations on the right hand side of Eq. (22) 
and Eq. (23) conserve energy, are accurate to order p ^, 
and will be useful when p ^ « 1 . These approximations 
could also have been derived directly from ray tracing. 
The extension of this type of analysis to three or more 
reflecting boundaries for all powers of p probably 
becomes impractically complex. On the other hand, ray 
tracing to generate all terms linear in p appears relative- 
ly straightforward. But, the result will be usefiil only 
when the interference effects are small compared to the 
effects of interest. 



4. Examples 

Reference [2] describes in detail an application of the 
equations given above; including a comparison of the 
derived results with experimental measurements of 
microwave heating of water in a microchannel. Figure 2 
shows a top and a side view of the coplanar-waveguide 
transmission line of [2]. It is comprised of metal lines on 
a borosilicate glass substrate and a rectangular poly- 
dimethylsiloxane (PDMS) cover attached to the trans- 
mission line by contact adhesion. Region 1 of the trans- 
mission line extends from the microwave probes on the 
left side of the transmission line to the left end of the 
microchannel. Region 2 extends from the left end to the 
right end of the microchannel. Region 3 extends from the 
right end of the microchannel to the microwave probes 
on the right side of the transmission line. Regions and 
1 extend from the microwave probes to the closest end of 
the cover, respectively. 

Figure 3 is a closer view of the part of Fig. 2 that 
shows the microchannel in the cover over the 
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Fig. 2. Top and cross-sectional views of a co-planar transmission line comprised of a cover over metal lines deposited on a substrate. 
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Fig. 3. Close up view of the removable cover of Fig. 2 showing that the cover contains a microchannel. 



coplanar waveguide. Figure 3 is not to scale, and the 
channel cross section, which is actually trapezoidal, is 
only shown schematically. Figure 4 plots the transmis- 
sion-line reflectance spectra measured with a vector 
network analyzer, as described in Ref. [2], both when 
the microchannel in Fig. 3 was empty, and when it was 
full of de-ionized (DI) water Figure 5 plots the trans- 
mittance spectra for the same configurations of the 
transmission line. Most of the features on the reflec- 
tance spectra are probably interference fringes caused 
by multiple reflections between the boundaries of the 
homogeneous regions in the transmission line. Any 
interference fringes present in the transmission data 
are much less pronounced than on the reflection data 
due to the different scale on the y-axes in these two 
figures. 

The metal lines in the transmission line were 500 nm 
thick, vacuum-deposited gold on thin chromium 
adhesion layers, Zy~z^~ 0.569 cm, and Zj = 0.362 cm. 



Over the frequency range of Figs. 4 and 5 borosilicate 
glass and PDMS are approximately lossless (do not 
absorb microwave power) and have dielectric constants 
not very different from air when compared with those 
of gold and water. Therefore, the first assumption 
needed to apply the above equations to the transmission 
line of Fig. 2 is that the transmission line properties 
in Regions 0, 1,3, and 4 in Fig. 2 are identical. This 
is illustrated by showing Region as a subset of 
Region 1 and Region 4 as a subset of Region 3 in 
Fig. 2. Thus, the three-region transmission line 
in Fig. 1 approximates the five-region transmission line 
in Fig. 2 for the analysis presented here. 

Figure 6 plots the empty- and full-channel absorption 
spectra of [2] 
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(29) 
(30) 
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Fig. 4. Reflectance spectra of the microfluidic transmission line 
described in [2]. 
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Fig. 5. Transmittance spectra of the microfluidic transmission line 
described in [2]. 
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Fig. 6. Absorption spectra of the microfluidic transmission line 
described in [2]. 



respectively, where subscript e refers to the empty- 
channel data, and subscript / refers to the full- 
channel data. Any interference fringes on the empty 
and full-channel absorption spectra below about 
25 GHz are very small compared to the differ- 



ence between the two spectra. This fact suggests that 
the approximate equations derived here, which ignore 
interference effects, can be used to calculate the frac- 
tion of the incident microwave power that is absorbed 
in the channel water with little error due to interference 
effects below this frequency. Features that may be 
interference fringes are more evident above 25 GHz, 
but they are still relatively small compared to the dif- 
ference between the absorption spectra, so the errors 
incurred will be relatively small, even in this spectral 
region. 

4.1 Difference in Absorbed-Power Fractions 

The simplest estimate of the fraction of the 
microwave power incident on the filled transmission 
line that is absorbed in the water in the microchannel is 
the difference between the fraction of the power 
absorbed when the channel is frill and when it is empty, 
which is given by 



AA = A^-A^ 



R.+T-R 



'/ 



■T 



(31) 



However, as will be shown below, this is only a 
lower bound to the fraction of the power absorbed in 
the microchannel water and can greatly underestimate 
the actual fraction absorbed. 

4.2 Reflection Only From the Left Boundary 
of Region 1 

To derive a more accurate approximation, further 
assumptions about the transmission line are required. 
For instance, it can be assumed that Eqs. (1-5) describe 
the empty device. In other words, Ri2 = R23=R3 = and 
p 5^ is the power reflection coefficient at the left 
boundary of Region 1 . Let R = R^ and T= T^ be the 
measured reflectance and transmittance for the empty 
transmission line in Eqs. (1) and (2). In this case, 
these equations can be combined to give the attenuation 
coefficient of the empty transmission line as 



a„ 



-1 



2[Zj + Zj + Z3 ] 



-In 



r„ "l 



1-^. 



(32) 



and A, 



and A.,„ can be calculated from 



Eqs. (3-5), and Eq. (32), with ai = a2= 0C2= a^. Here 
the numeric subscripts designate the region, sub- 
script e designates the empty transmission line, and 
subscript m designates absorption due to the metal con- 
ductors. 

Now let R=Rf and T=Tf be the measured 
reflectance and transmittance calculated from the 
S-parameters when the channel is filled with DI water 
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The next assumption required is that the transmission- 
line attenuation coefficient due to the metal conductors, 
does not change when the microchannel in the trans- 
mission line is filled with water This is a good approx- 
imation, provided that the addition of the water does 
not substantially change the field distribution on the 
surface of the metal conductors, even if it changes the 
distribution significantly in the middle of the channel. 
Regions 1 and 3 are still identical, and «! = aj = a„, but 
0^2 = cc„ + a„ where a„ is the portion of the attenuation 
coefficient of the transmission line in Region 2 due to 
the additional absorption by the water in the channel. In 
this case, Eqs. (1) and (2) can be combined to give 



«,., = - 



-1 



2z, 



In 



( t\ 1 (33) 

2aJZi+Z2+Z3]^ 



\-R. 



and Ayf^, A^^, Ajfi^, and A^f^ can be calculated from 
Eqs. (3-5), Eq. (32), and Eq. (33) with «! = ttj = a^ and 
0^2 = cc„ + a„. Again, the numeric subscripts designate 
the region, but subscript / designates the filled 
transmission line, and subscript w designates absorp- 
tion due to the water in the microchannel. 

Table 1 presents the results of applying Eqs. (32) and 
(33) to S-parameter data measured at 5.02 GHz [1]. The 
first and second sections of the table compare the 
measured and calculated parameters for the empty 
transmission line and the transmission line filled with 
DI water, respectively. First, notice that the calculation 
conserves power in that the sum of the reflected, trans- 
mitted, and absorbed fractions of the incident power is 
unity. Also note that R_f^^ R^ shows at least one para- 
meter changed when the transmission line was filled 
with water, and explains whyA^^f^^ Ay in Table 1 . 

Also in Table 1, notice that T^^T^ and ^2f'^^2e- 
A small fraction of this difference is attributable to 
violations of the assumptions on which the analysis 
was based. The main contributor is the large 
difference between the attenuation coefficients 
ttj = oc^ for the empty channel case, and 02 = a^ + a„ 
when the channel is filled with water In other words, 
filling the microchannel with water has doubled the 
absorption (attenuation) in Region 2, with about half of 
the power being absorbed in the metal conductors and 
the other half being absorbed in the water in that region. 

More precisely, the fraction ^2/m of the power 
absorbed in the Region 2 metal of the full-channel 
transmission line is only about 93 % as large as the 



Table 1. Measured and calculated parameters for the device 
described in [1] at 5.02 GHz. The subscript numbers correspond to 
the different regions. The subscripts e and/correspond to an empty 
channel and a channel filled with water, respectively, and the sub- 
scripts m and w correspond to the metal conductors and channel 
water 
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Parameter 




Value 


Unit 


Empty 




^11. 




-15.9956 


dB 
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^2ie 

Re 

Te 




-3.24646 
0.025144 
0.473537 


dB 






«„, 




0.240687 


cm"' 






^l™ 




0.233569 








^2™ 




0.118546 








^3™ 




0.149204 






«.+r,+^i™+^2„+^3™ 


1.000000 




Filled 




^W 




-22.6425 


dB 


Channel 




Rf 

Tf 




-3.97558 
0.005442 
0.400352 


dB 






a„ 




0.259523 


cm"' 






^ijh, 




0.238289 








^2Jh, 




0.110560 








Ay^ 




0.119212 








A,„ 




0.126145 






Rj+T/-A 


l>"'"^2/m"'"^2>"'"^3>i 


1.000000 





fraction Aj^^ absorbed in the empty-channel transmis- 
sion line. This difference is not the result of experimen- 
tal error or a violation of the assumptions used in deriv- 
ing the model. Instead, this difference is to be expect- 
ed. When the water is present, the microwave power 
decreases more rapidly with distance along the channel 
than when the channel is empty due to the addition of 
water absorption to metal absorption. Therefore, there 
is less power available for the metal in Region 2 to 
absorb when the channel is full than when it is empty. 
In other words, water absorption competes with metal 
absorption in Region 2 when the channel is full of 
water 

Figure 7 compares the amplitude-attenuation coeffi- 
cients for the metal conductors a„ and the channel 
water a„ that were calculated from the measured 
reflectance and transmittance spectra as illustrated in 
Table 1 . Inspection of Fig. 7 reveals that interference 
fringes are evident in both spectra, but are most obvi- 
ous in the a„ spectrum. These are almost certainly 
errors caused by ignoring interference among multiple 
reflections in the approximate equations. But, these 
errors are small compared to the major trends in the 
spectra. 
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Fig. 7. Microwave attenuration coefficient data for the metal 
conductors and the channel water for the transmission line described 
in [2]. 



4.3 Alternative Assumptions and Comparison of 
Results 

The model given in Eq. (31) of Sec. 4.1 can be 
rewritten as 



Af.=^ 



(34) 



in terms of the definitions given in connection with 
Table 1, where AA is defined in Eq. (31). Recall that 
Eq. (31) was derived by ignoring the competition 
between absorption by Region 2 water and metal. 
For comparison with other models considered in this 
paper, the model of Sec. 4.1 will be referred to as 
Model (a). Similarly, the model discussed in Sec. 4.2 
based on the assumption that only the reflection coeffi- 
cient at the left probe is non-zero, will be referred to as 
Model (c). As alternatives to these models, the reflec- 
tion coefficients might be assumed to be 

(b) zero because they are small, 

(d) non-zero only at the left end of the micro- 
channel, 



These possibilities will be considered next. For case 
(b), which assumes that all reflection coefficients 
are zero, the same equations used in Sec. 4.2 can be 
used with the assumption that R^ = Rf= 0. For case (h), 
Eqs. (1) and (12) can be combined to eliminate TiTjTj in 
these equations with the result 



2R + r 



-VF+4M^ 



(35) 



2R 



which can be substituted first into Eq. (32) to determine 
a„ and then along with a„ into Eq. (33) to determine 

Closed form solutions are not possible for the 
remaining equations, so they must he solved numerical- 
ly. For case (d), which assumes non-zero reflection 
only at the left end of the microchannel, the approxima- 
tions on the right hand side of Eqs. (7) and (1) can be 
combined to eliminate p to give. 



R T 

+ 



= 1, 



(36) 



which can then be solved numerically, first, for a^ in 
terms of R^ and T^ and then for a^ in terms of a„, R„, 
and T„. 

For case (g), which assumes non-zero reflection only 
at the right end of the microchannel, the approxima- 
tions on the right hand side of Eqs. (9) and (1) can be 
combined to eliminate p to give, 



R 

2 2 



= 1. 



(37) 



which can then be solved numerically, first, for a^ in 
terms of R^ and T^ and then for a„ in terms of a„, R„, 
and T„. 

For case (e), which assumes non-zero reflection only 
at both probes, the approximations on the right hand 
side of Eqs. (16) and (17) can be combined to eliminate 
p to give. 



2R 

1 + t' 



= 1. 



(38) 



(e) non-zero only at both probes, 

(f) non-zero only at both ends of the microchannel, 

(g) non-zero only at the right end of the micro- 
channel, 

(h) non-zero only at the right probe. 



which can then be solved numerically, first, for a^ in 
terms of R^ and T^ and then for a„ in terms of a„, R„, 
and T„. 

For case (f), which assumes non-zero reflection only 
at both ends of the microchannel, the approximations 
on the right hand side of Eqs. (24) and (25) can be 
combined to eliminate p to give. 
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2R 



Till + Tj] TiTjT, 



= 1, 



(39) 



which can then be solved numerically, first for a^ in 
terms of R^ and T^ and then for a„ in terms of a^, R„, 
and T„. 

Table 2 lists values of R^ and T^ obtained at 5 GHz, 
15 GHz, and 23 GHz [1,2]. Table 3 compares the 
results of applying the different models derived above 
to the data in Table 2. The approximation Ay^ = AA 
of Eqs. (31) and (34), which is the simplest estimate of 
the fraction of the microwave power absorbed in the 
water, is also shown in the table for comparison. Notice 
that this estimate is significantly smaller at all 
three frequencies than that based on the other approxi- 



mations discussed in this paper This is a general result 
because AA completely ignores the competition 
between the water and metal absorption processes, 
which reduces the total amount of power absorbed 
in the metal when the water is present, as discussed 
previously. 

Table 2. Measured transmittance-line reflectance and transmittance 
at three fi'equencies 



Frequency 



R 



^f 



V 



5 GHz 


0.025144 


0.473537 


15 GHz 


0.006318 


0.308259 


23 GHz 


0.028247 


0.244598 



0.005442 0.400352 
0.004613 0.203731 
0.028454 0.142578 



Table 3. Comparison of results fi'om the different models discussed. The results in the column headed by p,^ are based on the assumption that 
the reflection coefficient at the left probe was non-zero. Similarly, the resuhs in the column headed by p„ are based on the assumption that the 
reflection coefficient at the right end of the channel is non-zero, etc. The other quantities are defined in the text. 

















Non-Zero Reflection Coefficients 








KA 


all zero 


P,P 


Pic 


Plp = Prp 


Plc = Prp 


Prp 


Prp 


Condition 


Parameter 


(a) 


(b) 


(c) 


(d) 


(e) 


® 


(g) 


(h) 


5 GHz 


Pe 




0.0000 


0.0251 


0.0429 


0.0202 


0.0249 


0.0590 


0.0924 


Empty 


Km 




0.2492 


0.2407 


0.2346 


0.2354 


0.2322 


0.2289 


0.2169 


Channel 


^2.» 




0.1243 


0.1185 


0.1145 


0.1186 


0.1183 


0.1236 


0.1190 


5 GHz 


Pf 




0.0000 


0.0054 


0.0093 


0.0047 


0.0063 


0.0194 


0.0318 


Full 


aw 




0.2319 


0.2595 


0.2789 


0.2759 


0.2850 


0.2888 


0.3211 


Channel 


^2/» 




0.1147 


0.1106 


0.1076 


0.1087 


0.1074 


0.1079 


0.1029 




^2/» 


0.0929 


0.1068 


0.1192 


0.1279 


0.1274 


0.1319 


0.1191 


0.1524 


15 GHz 


Pe 




0.0000 


0.0063 


0.0152 


0.0058 


0.0097 


0.0263 


0.0589 


Empty 


a™ 




0.3928 


0.3902 


0.3871 


0.3884 


0.3858 


0.3834 


0.3720 


Channel 


^2.» 




0.1582 


0.1569 


0.1549 


0.1569 


0.1568 


0.1598 


0.1576 


15 GHlz 


Pf 




0.0000 


0.0046 


0.0111 


0.0044 


0.0089 


0.0426 


0.0917 


Full 


«». 




0.5720 


0.5744 


0.5778 


0.5757 


0.5742 


0.5488 


0.5230 


Channel 


^2/» 




0.1309 


0.1298 


0.1283 


0.1296 


0.1292 


0.1333 


0.1325 




^2/» 


0.1062 


0.1907 


0.1911 


0.1916 


0.1922 


0.1923 


0.1908 


0.1862 


23 GHz 


Pe 




0.0000 


0.0282 


0.0773 


0.0265 


0.0495 


0.1354 


0.2591 


Empty 


«» 




0.4694 


0.4598 


0.4425 


0.4512 


0.4346 


0.4209 


0.3694 


Channel 


^2.». 




0.1689 


0.1631 


0.1529 


0.1635 


0.1624 


0.1789 


0.1673 


23 GHz 


Pf 




0.0000 


0.0285 


0.0779 


0.0278 


0.0641 


0.2791 


0.4396 


Full 


«». 




0.7455 


0.7452 


0.7446 


0.7416 


0.7001 


0.4944 


0.3599 


Channel 


^2/» 




0.1325 


0.1279 


0.1198 


0.1278 


0.1262 


0.1579 


0.1517 




^2/» 


0.1018 


0.2104 


0.2072 


0.2015 


0.2101 


0.2032 


0.1855 


0.1478 
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The data at different frequencies presented in Table 3 
illustrate different relations among the empty- and full- 
channel transmission-line reflectances R^ and Rj^ The 
values of R^ and Rfars given in Table 2, but they are 
also given by p^ and Py^under (c) in Table 3, because this 
is the case where the only reflection occurs at the left 
probe. Relatively speaking, at 5 GHz, R^ is large 
(« 0.025) while iJ^is small (« 0.005) and ^2/w is small 
(< 0.15). At 23 GHz, R^ and R^ are large, R^ -R, is 
small, and A^f^ is large (> 0.15). At 15 GHz, both R^ 
and Rf and their difference are small, but ^j/w is large. 

Figure 8 compares the fraction ^j/w of the incident 
power that is absorbed in the microchannel water as 
calculated with the different models. Models (a) and (h) 
are plotted with squares, and the other models are 
plotted with times signs. The data points for (a) and (h) 
are connected with dashed lines, and the data points for 
(c) are connected with solid lines. An ellipse encircles 
a cluster of data points at 5 GHz, and another ellipse 
encircles a cluster of data points at 23 GHz. The cluster 
is so tight at 15 GHz that no ellipse is needed. The 
models that have at least one non-zero reflection coef- 
ficient on the left side of the microchannel (c - f) agree 
among themselves, significantly better than the other 
models. Also (b), which is based on assuming all reflec- 
tion coefficients to be zero, agrees much better with 
these than do the models that assume non-zero reflec- 
tion coefficients only to the right side of the micro- 
channel. The reasons for and implications of these 
facts, are discussed in the next section. 



5. Discussion and Conclusion 

Equations that describe the absorption of microwave 
power in the lossy components of a complex transmis- 
sion line have been derived for a number of simplified 
models of the general three-region transmission line. 
Examples of the use of these equations to estimate the 
fraction of the power absorbed in the water in a 
microchannel embedded in a transmission line were 
described. These examples required three assumptions 
common to all of the examples. First, it was assumed 
that each transmission line could be divided into three 
homogeneous regions extending from the probes on 
one side of the device to the probes on the other side of 
the device. Second, it was assumed that the metal con- 
ductors and channel water are the only lossy compo- 
nents of the transmission line. Third, it was assumed 
that the portion of the attenuation coefficient of the 
transmission line that describes absorption of 
microwave power in the metal conductors did not 
change when the channel was filled with water This is 
a good assumption if the presence of the water does 
not significantly change the field distribution on the 
surface of the metal conductors. 

For some of the examples, it was further necessary to 
assume that the reflections at three of the four region 
boundaries in the transmission line were zero, even 
though what appeared to be interference fringes on the 
transmission line reflectance spectra suggest that was 
not the case. The results derived for these examples are 
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Fig. 8. Comparison of the fraction of the incident microwave power absorbed by the 
water in the channel in the transmission line described in [2]. Models (a) - (h) are 
described in tbe text. 
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rigorous. For two of the examples, (e) and (f), it was 
assumed that reflections at two of the four discontinu- 
ities were zero. These models are the most realistic 
because they assume identical reflection coefficients at 
pairs of boundaries that are mirror images of each other 
and should therefore be identical in theory, and similar 
in practice. On the other hand, the results based on 
these models are accurate only to the extent that inter- 
ference effects are negligible, and then only to order p^, 
which will be more than adequate in many applications. 

At the three frequencies for which data are reported 
in Fig. 8, the results of four of the models (c - f) form 
tight clusters. Ellipses encircle these clusters at 5 GHz 
and 23 GHz, while the data in the cluster at 15 GHz 
cannot be resolved at the scale of the figure. Model (b) 
gives results that agree much better with the data in the 
clusters than do (a) and (h). It is reassuring in this con- 
nection that (a) and (h) are the least realistic models of 
the transmission line shown in Fig. 2. 

It has already been mentioned that (a) ignores the 
effect of the competition between the absorption by 
water and the metal conductors, and it only provides a 
lower bound to the fraction of the power absorbed in 
the water For the transmission line analyzed here, (h) 
requires an unrealistically large reflection coefficient at 
the right-hand probe to compensate for the absorption 
that is associated with two complete passes along the 
transmission line. To a lesser extent, (g) also suffers 
from this problem. 

The fact that five of the eight models compared here 
give substantially identical results for the fraction of the 
incident microwave radiation that is absorbed in a 
water filled channel under a range of reflection condi- 
tions suggests that these models can be used with some 
confidence. Specifically, they should provide reason- 
ably accurate results when the measured reflectance as 
well as the interference fringes on the measured 
reflectance and transmittance spectra, are small com- 
pared to the absorbed fraction of interest. The reason 
for this is, that the fraction of the microwave power 
absorbed in a transmission-line region as calculated 
from these models is not very sensitive to the detailed 
assumptions about which boundaries are reflecting. 
On the other hand, the calculated values of the re- 
flection coefficients are quite sensitive to the 



same assumptions. Nevertheless, if it is not obvious 
from physical reasoning or examination of experimen- 
tal reflection and transmittance spectra which of these 
reflection models best describes a particular complex 
transmission line, then it is probably a good idea to 
compare the results of a number of models, as was done 
here. 

In the general case, three-region devices will have 
four non-zero, non-identical reflection coefficients. 
Similar approximate equations that are accurate when 
the products of all combinations of two of the power 
reflection coefficients are small, can be easily derived 
for this case, but more independent measurement 
results will be needed to use these equations as illustrat- 
ed here. The results could be obtained by performing 
the S-parameter measurements with the same PDMS 
cover located in different positions along the axis of the 
coplanar waveguide to change the dimensions of Zj and 
Z3. This will change the fractions of the power returned 
to Boundary A by the reflections from Boundaries B, C, 
and D, and will change them in different ways for each 
boundary. 
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